Force-extension relation of cross-linked anisotropic polymer networks 
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Cross-linked polymer networks with orientational order constitute a wide class of soft materials 
and are relevant to biological systems (e.g., F-actin bundles). We analytically study the nonlinear 
force-extension relation of an array of parallel-aligned, strongly stretched semiflexible polymers with 
random cross-links. In the strong stretching limit, the effect of the cross-links is purely entropic, 
independent of the bending rigidity of the chains. Cross-links enhance the differential stretching 
stiffness of the bundle. For hard cross-links, the cross-link contribution to the force-extension relation 
scales inversely proportional to the force. Its dependence on the cross-link density, close to the 
gelation transition, is the same as that of the shear modulus. The qualitative behavior is captured 
by a toy model of two chains with a single cross-link in the middle. 

PACS numbers: 82.70.Gg, 64.70.Md, 61.43.Fs, 87.16.Ka 



I. INTRODUCTION 

Anisotropic networks are abundant in nature and 
among man-made materials. For example, nematic 
elastomers PQ or actin- myosin assemblies in cells [2]. 
These networks are constructed from anisotropic build- 
ing blocks, nematogens such as stiff rods or semiflexible 
polymers which are well modelled by wormlike chains 
(WLCs) [3] . Isotropic networks of WLCs have been stud- 
ied extensively in experiment [4i, driven by the interest in 
biological networks such as the cytoskeleton or the extra- 
cellular matrix. Theoretical approaches have dealt with 
entangled solution [5] as well as chemically cross-linked 
networks of WLCs [SJ [7] . The former build on the tube 
model whereas the latter generalise concepts from rub- 
ber elasticity. Aligning linkers can give rise to a variety 
of morphologies which have been studied by means of a 
generalised Onsager approach [5] as well as within a mi- 
croscopic model OUS]. Of particular interest are bundles 
of filaments which occur in a broad range of cytoscele- 
tal structures and show much richer elastic behaviour 
than the usual WLC Q3JCE31H1]. Tw0 parallel-aligned, 
stretched semiflexible filaments cross-linked by a motor 
cluster have been used as the minimal elastic element of 
an active gel [15] , 

Here, we consider an anisotropic network of WLCs 
which have been aligned along a preferred axis chosen 
as the z-direction. The alignment is not due to cross- 
links which we model as springs. Instead possible mech- 
anisms for alignment are a nematic environment, pulling 
forces, grafting surfaces or the Onsager mechanism. A 
sketch of such a network is shown in Fig. [T] If the align- 
ment is strong, the WLC model can be replaced by a 
weakly bending chain (WBC) as first suggested by Marco 
and Siggia [16] for DNA molecules. The advantage for 
an analytical approach is enormous because the single 
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FIG. 1: network of aligned wormlike chains under pulling 
force 



chain model is Gaussian. Our focus in this paper lies 
on the force-extension curve of a randomly cross-linked 
anisotropic network of strongly aligned filaments. We 
first consider a toy model, consisting of two cross-linked 
filaments. The model allows us to disentangle the contri- 
butions to the effective extension which are due to either 
bending stiffnes or cross-links. Subsequently, we analyse 
the effects of cross-links for a macroscopic network. 
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II. MODEL 



Our starting point is the energy of a stretched WLC 



in terms of the tangent vector t(s) 



_ rfr(s) _ i dr± d 
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H[t(s)] 
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Here k denotes the bending stiffnes which is related to 
the persistence length L p via L p = 2n/{{d — \)ksT), 
where d is the dimensionality of the embedding space. 
The pulling force is denoted by /, and < s < L is 
the arclength. The local inextensibility constraint of the 
WLC is expressed by the condition |t(s)| = 1. We assume 
that the chain is strongly stretched so that tilting of the 
tangent vector away from the z-axis is small and we can 
use the approximation 




FIG. 2: Two aligned chains with one cross-link in the middle 



should be 



1 /dr ± 



2\ ds 



(2) 



leading to the weakly bending model introduced by 
Marco and Siggia [16] 



1 - d 2 r ± (s) y f f L d {dr x (8)\* 
ds? ) + 2 J 

(3) 

The central quantity of interest is the extension of the 
chain under an applied force /, which in the weakly bend- 
ing approximation is computed from the thermal fluctu- 
ations transverse to the aligning direction: 

(z(L)-z(O)) = fds(f s ) 



ds 



2 \ ds J 



(4) 



III. TOY MODEL: 2 CROSS-LINKED CHAINS 

Before addressing the full problem of a randomly cross- 
linked array of aligned chains, we discuss the much sim- 
pler case of two strongly stretched chains in two dimen- 
sions with one cross-link in the middle, see Fig. [2] In two 
dimensions, rj^(s) = y(s) and the Hamiltonian reads 



H 



Wobi(s)]+Ko[yi(s)] + |( yi (i/2)-2/2(i/2)) 2 (5) 



The cross-link is modeled as a harmonic spring of stiff- 
ness g. For simplicity we impose hinged-hinged boundary 
conditions: yi(0) = yi(L) = 0, 2/2(0) = 2/2 (£) = D and 
y'{(0) - y"(L) = 0, y 2 '(0) = y$(L) = for the two 
chains which are a distance D apart (the prime denotes 
derivative with respect to s) . According to the boundary 
conditions that we use, the eigenfunction representation 



1=1 
00 

y 2 (s) - D = ^2 B l sin(<? ; s) 
1=1 

and wavenumbers are restricted to values qi — j} , I G Z. 
These eigenfunctions diagonalise the Hamiltonian of the 
weakly bending chain, whereas the cross-link gives rise 
to a term which is quadratic in the amplitudes but not 
diagonal. Introducing vectors 

r = (A 1 ,B 1 ,A 2 ,...) , 

u = g 1/2 {sm( qi L/2), - sm( qi L/2), sm(q 2 L/2), ...) 



and matrices 



C = 



L 



fq 2 i^qf + f) 

q\{nql + f) 

q^qf + f) 



V 



J 



the Hamiltonian is rewritten as 



(6) 



l,m=l 1=1 

with G;. ln = Ci_ m +uiu m . The matrix G is easily inverted: 



G" 1 = C" 1 - 



C- x uu T C- x 
l + u T C- l u ' 



(7) 



so that we can compute the force-extension curve, say of 
chain 1, exactly 



L 



(zi(L)) = L--J2qf(A 2 ) 



(8) 



Our model has three characteristic energies: the bend- 
ing energy, k/L, the work done by the external force, fL 
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and the thermal energy, T (we set ks = 1)- Use of the 
weakly bending approximation requires k/L»T (which 
is equivalent to L p /L S> 1) or fL » T. This leaves us 
with one free parameter, x := /L/(re/L), namely the ra- 
tio of work done by the external force to bending energy. 
This dimensionless quantity can also be interpreted as 
the squared ratio of two lengthscales: the total contour 
length L to the length yj re// over which the boundary 
conditions penetrate into the bulk (i.e., the size of a link 
in an effective freely-jointed chain) [17]. Actually, we 
have two additional lengthscales, the distance between 
chains, D, which we assume to be negligible and the 
length of the cross- link (cx 1/y/g) or alternatively the 
energy of a cross-link relative to the work done by the 
pulling force, gL/ f. 

The result for general x and cross-link strength g 
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XL 



(9) 



can be decomposed into a contribution, Am St which is 
characteristic for a weakly bending chain and well-known 
from the work of Marko and Siggia |16j and a contribu- 
tion due to the cross-link, Axl 



A 



XL 



u (x) 
u 2 (x) 



gT u 2 (x) 
16/ 2 u (x) 



x- 2tanh( v / i/2), 
cosh( v / x/2)) 2 - 6 tanh(y/x/2) 



In the limit of large x, the cross-link contribution is given 
by 



A XL = 



8/ 2 ^ 2/ 



(10) 



For soft cross-links Axl falls off as /~ 2 , whereas in the 
limit of hard cross- links, gL/f — > oo, we find Axl = 
— T/(4/L). In any case, the contribution of the cross- 
link is subdominant in the limit of strong pulling force. 
Note that Axl is independent of re, which shows us that 
the contribution from the cross-link is purely entropic. 
This contribution could in fact be evaluated from a di- 
rected polymer model which does not involve any bending 
rigidity. 

Collecting the leading terms for hard cross-links and 
strong pulling force, yields 



L 



- 1 = - 



T T 
(4re/)V2 _ ifL 



8L?' 



(11) 



Here we have restored the term due to a finite distance, 
D, between the chains, which gives rise to a geometric 
reduction in length due to the cross-link. This term is 
presumably unimportant in a network, where D/L is ex- 
pected to be small. This term will be neglected in the 
following. 

The weakly bending approximation is satisfied as long 
as the r.h.s. of Eq. ^ is small. We point out that, in the 
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FIG. 3: Force-extension curve for L/L p — 1; red (black) curve 
without (with) cross-link 
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FIG. 4: Inverse differential stiffnes for L/L p 
curve without (with) cross-link 



1; red (black) 



case of L p ^> L, this approximation is fulfilled even with- 
out having very large x (strong stretching). We show the 
relative extension as a function of pulling force in Fig. [3] 
for hard cross-links. As discussed, cross-linking enhances 
the extension due to the reduction of thermal fluctua- 
tions, but the effect becomes less and less pronounced in 
the strong stretching limit. It is also of interest to con- 
sider a situation, where a strong pulling force has been 
applied and subsequently the change in extension in re- 
sponse to a small change in the pulling force is measured. 
This response is determined by the differential stiffness 
df /d{z\). In Fig. [4] we show the inverse stiffness for the 
same set of parameters as in Fig. [3] 

As expected the cross-link enhances the differential 
stiffness. If the chains are already strongly stretched, the 
cross-link has little effect. However for weakly stretched 
chains the enhancement is considerable. 



IV. RANDOMLY CROSS-LINKED NETWORK 

In this Section, we are going to compute the force- 
extension relation of a randomly cross-linked ensemble 
of oriented chains approximately. We are guided by the 
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result for the cross-linked pair of chains, which consists 
of a single chain contribution, ~ T '/ '(nf) 1 ^ 2 , and a con- 
tribution due to the cross-link, ~ T/(fL), which is in- 
dependent of the bending rigidity k. We decompose the 
calculation for the network accordingly: In the uncross- 
linked network, the force-extension relation is determined 
by the single chain contribution. To assess the effect of 
cross-linking, we compute the free energy difference of 
the cross-linked network relative to the uncross-linked 
system. 

Our starting point is the Hamiltonian 



u 



N 

£ 

i=l 



U [vi(s)] + H cv + Hxrmk(C M ) 



(12) 



with H given in Eq. ([3]). To simplify the notation, we 
have dropped the subscript ij_(s) — > Ti(s), such that 
denotes the transverse excursion of chain i. The excluded 
volume interaction 



^v = E? /^(r^)-!^)) (13) 

is introduced to balance the attractive interactions due 
to cross-linking. The latter are modelled by harmonic 
springs 



M 



(14) 



e=l 



where Cm '■= {*e>,7e> s e} is a quenched configuration of 
M cross-links connecting polymers i e ,j e at arclength s e . 

The partition function of the cross-linked system rel- 
ative to the uncross-linked melt for a specific realization 
of cross-links, Cm, reads 



Z(C M ) = (exp - 



'Hxi inkjCfcf) 

T 



(15) 



Here (...) denotes the thermal average over all polymer 
configurations with the Boltzmann weight exp(— %' /T) 

of the uncross-linked melt H' = 2i=i ^o[ r i(s)] + W ev . 

Physical observables can be calculated from the 
quenched-disorder averaged free energy, AF = — T\ha.Z], 
where [...] denotes average over all realizations of random 
cross-links. We assume that the number of cross-links 
can vary and a realization with M cross-links follows the 
Deam-Edwards distribution 1181: 



P(C M ) cx 



( M 2 A 



1 



M\ \2N(2ira 2 ) 



Z{C M ) , 



(16) 



where a 2 = T/g. The parameter /i 2 controls the aver- 
age number of cross-links per polymer and the physical 
meaning of this distribution is that polymer segments 
close to each other in the melt have a high probabiblity 
to be linked. Of particular interest is the derivative of 



the free energy change due to cross-linking with respect 
to the pulling force 



dAF 



(17) 



which yields the mean extension per chain i ^2i[(zi(L) — 
Zi(0))] relative to the uncrosslinked melt. 

The above model is expected to have a gelation tran- 
sition P35] at a critical cross- link concentration /i 2 ~ 1. 
As far as the force-extension relation is concerned, we 
expect to find the single polymer contribution (Marko- 
Siggia) below the gelation transition and a correction due 
to cross-links above it. The latter will be computed from 
the directed polymer model (k = 0), because it is dom- 
inated by the long wavelength transverse excursions of 
the polymers which are correctly captured by the second 
term in Eq. ^ . A similar mechanism underlies the well- 
known observation |20j that the tranverse fluctuations of 
a strongly stretched wormlike chain are independent of 

AC. 

We compute the free energy difference, AF, for a net- 
work of cross-linked directed polymers, confined between 
two planes with their endpoints free to slide on them. 
This calculation is analogous to our previous work on di- 
rected polymers |19) and some details are layed out in the 
Appendix. We point out that our calculation is restricted 
to the vicinity of the gel point. Adding the single chain 
contribution and denoting the distance from the gelation 
point by e = [i 2 — 1, we find for the total mean force- 
extension 



Lj2[{ Zi (L)- Zi (0))]-l = 



NL 
Ams = 



T 



(4«/)Va (2fL) 
e 3 T 
3 fL + 3f 2 /g- 



(18) 



The above result is quite remarkable in several re- 
spects. Whereas the contribution due to bending shows 
the l/v7 behaviour typical for WLC, the cross- link con- 
tribution is proportional to 1// (hard cross- links) , which 
is characteristic of freely jointed chains. For strong 
stretching, i> 1, the cross-link contribution has quali- 
tatively the same dependence on the pulling force, /, and 
the cross-link strength, g, as the corresponding contribu- 
tion in the two-chain toy model of the previous Section. 
There is a singular contribution to the stretching stiff- 
ness at the gelation transition, which has the same scal- 
ing as the shear modulus, namely ~ e 3 . In analogy to the 
behaviour of the shear modulus in the well cross-linked 
regime [5T] , we expect the contribution of the cross- links 
to the force extension relation to scale as the density of 
cross- links, i.e. to be of the form A XL ~ ^ 2 T/(fL). The 
excluded-volume interaction, which is included in the 
gelation theory in order to prevent collapse of the system 
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upon cross-linking, does not affect the force-extension 
relation in the strong stretching regime. Replacing the 
excluded- volume interaction of each chain with its neigh- 
bours by an effective harmonic "cage" [35], one notices 
that for strong pulling forces the effect of the "cage" be- 
comes negligible [23] . 



V. CONCLUSIONS-OUTLOOK 

In conclusion, we have calculated the effect of random 
cross-links on the force-extension relation of an array of 
parallel-aligned WLCs. Our calculation is restricted to 
the strong stretching regime, close to the gelation transi- 
tion. The main result is a contribution to the nonlinear 
force-extension relation which scales as ~ 1// for hard 
cross-links and suppresses the thermal flucuations thus 
stiffening the (thermal) stretching modulus. Hard cross- 
links are sufficient in the large L limit as the cross-link 
size a enters through TL/(a 2 f). Our result is based on 
a replica field theory of randomly cross-linked directed 
polymers originally developed in Ref. |19j . Remarkably, 
apart from numerical prefactors, the effect of cross-links 
is the same as in the case of a simple two-dimensional 
model with two chains and a single cross-link in the mid- 
dle. 

An interesting extension of this work would consider 
cross-links which are non-local in the z-direction. That 
would imply discontinuities in the tension of the involved 
chains as discussed in Ref. [15] . This problem will be 
addressed in a future publication. 



(ro, ...,r„). With that, we can express 



[Z n ] = %±i with 



Z n+1 = ( exp 



( M 2 



\2NL(j) J 



Jo 




In the effective replica partition function Z n +i, the av- 
eraging (•••)„+! w is done with the statistical weights 
"Ho + T~Lcv, and convenviently, in this form, the dis- 
order average has not to be taken into account any- 
more. The degrees of freedom are the replicated parti- 
cle positions {ri{z)}i=i,...,N and the definition A(x) := 
cxp(— x 2 /(2a 2 )) resembles the interaction of the cross- 
linkers. 

In this form the Hubbard-Stratonovich transformation 
can be used to change the degrees of freedom from the 
particle positions to a (replicated and Fourier-space) den- 
sity field (l(9,s) = (sEi^ (>) >: 



-"n+1 



(22) 



Here, the effective replica Hamiltonian /„ + i(Sl) is given 
by: 



/«+i(fi) = /o + 



6> 2 
~2L~ 



V A(g) f L &s\n{q,s)\ 2 -\ni 

;r nDC JO 

(23) 



ijGHRS 



with the Fourier transform A(q) of A(i) and a single 
polymer partition function: 
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VI. APPENDIX 



Here, we outline the calculation of the free energy 
AF = -T[kiZ(C M )], Eqs. (p) and (pi of a disorder 



averaged randomly cross-linked network relative to the 
uncrosslinked state. For that, we start with the replica 
trick: 



AF 
~T~ 



= [In Z] = lim 



[Z n ] - 1 



(19) 



As usual [IH] this leads to replicated, D(n + 1)- 
dimensional vectors denoted by a hat, e.g. f — 



3= fvf(s)exp(^f- A (<?) fdsW,')^ 

J \ geHRS ^° 

(24) 

HRS stands for the higher replica sector, the set of q- 
vectors with at least two non-zero replica components. 
As been done several times before [TH [TS] , the excluded 
volume interaction is assumed to be strong enough to 
make the network incompressible. Hence density fluc- 
tuations, which would be represented by fl(q, s) with q 
having a non-zero component in only one replica, do not 
appear in the Hamiltonian / n+ i(f2). /o in ( |23[ ) is an 
unimportant contribution, which does not depend on f2. 

As next step, we perform the saddle point approxima- 
tion of ( 22 ) . The saddle point value of SI is given by: 



ft(q, s) = <2<S qo+ ... +q „, 



dev(e, 



s ) exp 



2 



(25) 



Here the gel fraction Q is the fraction of chains which are 
localized, which means they cannot traverse the whole 



G 



sample, but perform fluctuations around a preferred po- 
sition. The localization lengths £ quantify the extent 
of these in-plane fluctuations, which can depend on the 
height s in the sample. Their probability distribution 
■p(£ 2 ,s) has been determined in [T§] . 



We now plug the saddle point value £1 into Eq. ( 22 ) and 



restrict ourselves to the vincinity of the gelation transi- 
tion, i. e. small gel fractions Q. Bearing in mind that 
il oc Q, we can expand Eq. (24 1 in powers of Q and 



easily perform the functional integral over f(s): 



3 = 1 + A 4 Q 



2! 



-3(2) 



3! 



3(3) +<W 



(26) 



Here the coefficients 3(2) and 3(3) are given by: 



ijGHRS ^° 

a2 



I? 



Q 2 



x exp 



2/ 



\S2 - Si 



(27) 



3(3) = Yl ^ 1+ g 2+ g 3 ,6 A ('?i) A (92)A(g3) 

L dSlds 2 ds 3 H(gl, Sl)^(g2, S2)^(g3, S 3 ) 
L 3 Q 3 

( '*p( ~7T7 (9l?2|si-S 2 | +gl93|si-S3| + g293|s2-S3| 







2/ 



(28) 



For simplification, we present the calculation for hard 
cross-links, i. e. a = 0. However, the extension to ar- 
bitrary a is straightforward. Also, while we present the 
calculation for a three-dimensional system, the general- 
ization to arbitrary dimension is possible. The result for 
these generalizations is shown at the end. 

We now perform the sums over q. Assuming that the 
surface area A of the sample in the in-plane directions 
is large compared to microscopic details of the network, 
these sums can be changed to integrals. We obtain up to 
linear order in n: 



3(2) =l-n / d(£?,«i)d(£f,« 2 ) 

'g+g + |*l-*2|// 



In 



.4 



(29) 



3 (3) =l-n I d(tf,«i)d(&s 2 )d(&« 3 ) 
I In 



an + &e 3 + e 2 e 3 \ + 2c . (1 + (30) 



A 2 



Here, e :— fi 2 — 1 is the distance from the sol-gel tran- 
sition. It is related to the gel fraction by Q = 2e + 
C(e 2 ). For a better readability, we defined J d(£, 2 ,s a ) :— 

\ ds a J d£ 2 V{t; 2 , s Q ), and c = 1 + ln(27r) is a numer- 
ical constant. 

In a similar fashion, the sum over q can be performed 
in the second term of Eq. (23 1. 



^ £ A(q) f ds\ms)\ 2 = l-n f ^ 



(31) 



As one can see from Eq. (19), (|20|) and (22), the dis- 



order averaged free energy density is - in saddle point 
approximation - the term of /„ + i(J7) linear in n: 



l^AF = df n+1 
NT dn 



(32) 



71=0 



Hence, using the results (29)-fl31b, we can now recom- 



pose the free energy up to third order in e. With spatial 
dimension D + 1 and cross-link length a, the we obtain: 



1 dAF 
DNT df 



L 



3 Lf + 3f 2 a 2 



+ 0(e 4 ) 



(33) 



[1] M. Warner and E. M. Terentjev, Liquid Crystal Elas- 
tomers, (Oxford University Press, New York, 2003). 
[2] G. M. Cooper, The Cell, (Sinauer, Sunderland (MA), 



2000). 

[3] N. Saito, K. Takahashi and Y. Yunoki, J. Phys. Soc. Jpn. 
22, 219 (1967). 



7 



[4] A. R. Bausch and K. Kroy, Nature Phys. 2, 231 (2006). 
[5] D. C. Morse, Macro-molecules 31, 7030 (1998); 31, 7044 

(1998); 32, 5934 (1999). 
[6] D. A. Head, A. J. Levine and F. MacKintosh, 91, 108102 

(2003). 

[7] J. Wilhelm and E. Frey, Phys. Rev. Lett. 91, 108103 
(2003). 

[8] I. Borukhov, R. Bruinsma, W. M. Gelbart, and A. J. Liu, 
Proc. Natl. Acad. Sci. USA 102, 3673 (2005). 

[9] P. Benetatos and A. Zippelius, Phys. Rev. Lett. 99, 
198301 (2007). 

[10] M. Kiemes, P. Benetatos and A. Zippelius, Phys. Rev. E 
83, 021905 (2011). 

[11] M. M. A. E. Claessens, M. Bathe, E. Frey, and A. Bausch, 
Nature Materials 5, 748 (2006). 

[12] P. M. Goldbart, H. E. Castillo and A. Zippelius, Ad- 
vances in Physics 45, 393 (1996). 

[13] O. Lieleg, M. M. A. E. Claessens, C. Heussinger, E. Frey, 
and A. Bausch, Phys. Rev. Lett. 99, 088102 (2007). 

[14] C. Heussinger, M. Bathe, and E. Frey, Phys. Rev. Lett 



99, 048101 (1998). 
[15] T. B. Liverpool, M. C. Marchetti, J.-F. Joanny, and J. 

Prost, Europhys. Lett. 85 18007 (2009). 
[16] J. F. Marco and E. D. Siggia, Macromolecules 28, 

8759 (1995). 

[17] P. Benetatos and E. M. Terentjev, Phys. Rev. E 81, 
031802 (2010). 

[18] R. T. Deam and S. F. Edwards, Proc. Trans, R. Soc. 
(London) A 280, 317 (1976). 

[19] S. Ulrich, A. Zippelius and P. Benetatos, Phys. Rev. E 
81, 021802 (2010). 

[20] P. G. de Gennes, in Polymer Liquid Crystals, edited by 
A. Ciferri, W. R. Kringbaum and R. B. Meyer (Academic 
Press, New York, 1982) Chapter 5. 

[21] S. Ulrich, X. Mao, P. M. Goldbart and A. Zippelius, Eu- 
rophys. Lett. 76, 677 (2006). 

[22] D. Ertas and D. Nelson, Physica C 272, 79 (1996). 

[23] J. Wang and H. Gao, J. Mater. Sci. 42, 8838 (2007). 



